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Abstract 

^D ■ We discuss a new single parameter quasi-particle model and study the thermody- 

namics of (2+l)-flavor quark gluon plasma (QGP). Our model with a single parameter 
explains remarkably well the lattice simulation results of Fodor et. al. [1] . 
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1 Introduction : 

Quasi-particle model (qQGP) of quark gluon plasma (QGP) was first proposed by Peshier 
et. al. [2] to explain the non-ideal equation of state (EoS), observed in lattice gauge theory 
simulations (LGT) [3]. At finite temperature, instead of real quarks and gluons with QCD 
(quantumchromodynamics) interactions we may as well consider the system to be made up 
of non-interacting quasi-particles with thermal masses, quasi-quarks and quasi-gluons, and 
study the thermodynamics. Quasi-particles are quanta of plasma collective modes excited by 
quarks and gluons through QCD interactions. Initial quasi-particle model was found to be 
thermodynamically inconsistent [4] and also not able to fit the more recent LGT results [5]. 
Gorenstein and Yang reformulated the statistical mechanics (SM) to solve the inconsistency, 
but end up with an extra undetermined, temperature dependent terms in the expressions 
for pressure and energy density, which need to be phenomenologically chosen. However, 
as we discuss here [6], above reformulation of SM is not needed and standard SM may be 
applied to qQGP without any extra phenomenological terms. There is no TD inconsistency 
in our new qQGP model. Peshier's model with reformulated SM by Gorenstein and Yang 
has been studied by various groups [7, 8, 9, 10, 11] with different expressions for thermal 
masses, effective degrees of freedom, so on. Note that all of above works are based on the 
reformulation of SM of Gorenstein and Yang, which in fact based on mathematical identities 
involving derivatives with respect to temperature and chemical potentials, used to redefine 
average energy density {e) and number density (n) respectively. As we have shown recently in 
Ref. [6], we may skip this TD inconsistency problem and instead use the original definition 
of e and n, and making use of TD relations we may get all TD quantities. As a specific 
example, here we discuss (2-|-l)-fiavor QGP, studied by Fodor et. al. [1] using LGT. 

Of course, there are other models like HTL (hard thermal loop) [12], recent FMR (fun- 
damental modular region) gas [13] etc. based on QCD perturbative and non-perturbative 
calculations, but fails to fit LGT results near to the transition temperature Tc. At the 
same time phenomenological models of QGP, based on plasma theory with QCD inputs like 
SCQGP (strongly coupled quark gluon plasma) [14], our present qQGP model seems to fit 



remarkably well the LGT results with minimum number of parameters. All other qQGP 
models, field theoretical models [15], Ploykov loop models [16] also fit LGT results, but by 
adjusting 3 or more parameters. Of course, we know that near T = Tc, region of phase 
transition or cross over, it is a low energy phenomena and hence QCD can not be solved by 
analytical methods like perturbation theory because the coupling constant ctg is not small 
enough. Probably we need to formulate phenomenological models, just like in the case of 
hadron spectroscopy, to study TD of QGP near Tc- 

2 Phenomenological Model with /i = 0: 

QGP at thermodynamic equilibrium consists of interacting quarks and gluons which exhibits 
collective behaviour. Our basic assumption is that this system may be replaced by a system 
of non-interacting quasi-particles with quantum numbers of quarks and gluons. These quasi- 
particles have additional thermal masses which are equal to plasma frequencies. Here we 



differ from other qQGP models where, for example, the thermal mass was taken to be a/3/2 
times the plasma frequency. A general expression for thermal mass or polarization tensor 
is very complicated expressions which is a function of momentum and frequency. Only at 
high momentum limit it approaches a simpler form which on further approximations reduces 
to above form. In view of such a drastic approximation and since we use phenomenological 
model one may as well take rrith ~ ujp- In fact with this relation, we get better result 



than with rrith ~ J^/'^^ujp- Further important point is that the above dispersion relation 
is obtained using perturbation methods with temperature dependent density distribution 
function appropriate to ideal system. Then one formulates TD of a non-ideal system. In 
principle, this must be carried out in a self-consistent manner as discussed in Ref. [17], 
where, for example, density expression is an integral equation since Up depends on density. 
So we need to solve an integral equation self-consistently to get the the density. Here we 
avoid all above complications and as a phenomenological input, we assume that rrith = ^p- 
Following the standard procedure of statistical mechanics [18], the grand partition func- 



tion is given by, 

s,r 

where the sum is over energy states Er and particle number states Ns- a and (3 are defined 
as a = —fi/T and /3 = 1/T. Next on further simphfications and taking thermodynamic 
hmit, we get, 

oo 

q = lnQG = Tj2H^^^^~^'') ' (2) 

A;=0 

where q is called q-potential and =F for bosons and fermions. z = e^/^ = e~" is called 
fugacity. e^ is the single particle energy, given by, 

where k is momentum and m is the total mass which contains both the rest mass and 
thermal mass {rrith) of particles, rrith may depend on temperature T and chemical potential 
/i depending on QGP system. All other qQGP models starts from q-potential to get pressure 
with correction term to avoid TD inconsistency and proceed to evaluate other TD quantities. 
However, we avoid this trap of inconsistency problem by adopting a different procedure. We 
know that, for grand canonical ensemble, TD quantities like energy density and number 
density may be obtained by taking ensemble averages. That is, the average energy U is 
given by. 



U =< Er >-- 
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Note that the partial differentiation with respect to f3 above is just a mathematical trick 
to express U in terms of sum over single particle energy levels, e^, making use of Eq. (2). 
While differentiating, indirect dependence oi f3 = 1/T in the fugacity, z, and thermal mass, 
mth(T,fi), must be ignored by definition. Otherwise, we will not get back < Er >■ Now, on 
taking continum limit and after some algebra, we get. 
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where gf is the degenarcy and equal to Qg = 16 for gluons and equal to 12 n/ for quarks, n/ 
is the number of flavors with same mass. Ki and K2 are modifled Bessel functions of order 
1 and 2 respectively. 

Let us now consider our main topic, (2+l)-flavor system, studied by Fodor et. al. [1], 
using our model. It is a QGP with two light {u) and one heavy (s) quarks along with 
gluons. Let us flrst consider the case with zero chemical potential and take z = 1. Hence 
we get the energy density, expressed in terms of e(T) = e/Ss-, for the quark gluon plasma of 
quasi-partons is 
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where Eg is the Stefan-Boltzman gas limit of QGP, which may be obtained by taking high 
temperature limits of Eq. (4) for gluons and quarks separately and adding them, rrig is 
the temperature dependent gluon mass {rrith), which is equal to the plasma frequency, i.e, 
{2Nc + Hf). All quarks have both the thermal mass as well as the rest mass 
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and hence, the total mass may be written as 
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following the idea used in other qQGP models for the system with flnite quark masses. Only 
the difference is that our rrith is equal to the plasma frequency due quarks alone. That is. 
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Uf. irLgQ is the rest mass of up or strange quark, g"^ in thermal masses is 
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related to the two- loop order running coupling constant, given by, 

m = Stt / _ 3(153 -19n;)ln(21n(T/Ar)) ' 

"'^ ^ (33 - 2n^) ln(T/AT) V (33 - 2n/)2 ln(T/Ar) 



(7) 



where At is a parameter related to QCD scale parameter. This choice of as(T) is motivated 
from lattice simulations. Using thermal masses with above as, we can evaluate the e(T) from 



Eq. (4). Note that the only temperature dependence in e(T) comes from as{T), which has 
the same form as that of lattice simulations [19] with At as a free parameter, n'^^-^ and nl-^-^ 
are effective number of quarks which differ from 2 and 1, respectively, because of finite masses 
of up and strange quarks. However, handling of finite rest masses of quarks are different 
in LGT studies of Fodor et. al. and Bielefeld group [19]. Bielefeld group carried out the 
simulation with the ratio rnqo/T equal to constant and it is straight forward to calculate rig-^^ , 
where as Fodor et. al. used constant rriqQ and it is not clear how to calculate Ug'^^ . Hence 
as done in Ref. [20] we take in our calculations n^^^ = 2, n^^^ = 0.5 and n^ = 2.5. Similar 
normalization need to be made to fit the LGT results in other models also, either multiplying 
LGT data with a factor 1.1 [9] or model's data with .9 [11], or sometime using massive gluon 
[21] so on. From e, we may obtain pressure P by using a TD relation e = T^ — P for /x = 
system and we get 

P Po fT ,^£(T) 
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where Pq and Tq are pressure and temperature at some reference points. Results are presented 
in Fig. 1 along with LGT results. Note that earlier this phenomenological new qQGP model 
with a single system dependent adjustable parameter explained very well the LGT results 
of Bielefeld group [19] on QGP system with massless quarks as discussed in Ref. [6]. 

Another TD quantity which may be obtained from ensemble averaging is the number 
density, given by, 

<N>= h-d:hl = -^^^QG = z^\nQa = Y: ' /^.. ' (9) 
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and using the definition n =< N > /V, where V is the volume. 

3 Model with finite /i : 

Let us next consider (2+l)-flavor system with finite fi and LGT results are available for 
quark density Uq or baryon density (n^) and difference in pressure from /x = case (AP = 
P{T,fi) — P{T,fi = 0)). Here n is the quark chemical potential which is one third of the 
baryon chemical potential. The present LGT results are with fis coming from up quarks 



only and hence we need to consider only up quark density. Using the expression for the 
number density, Eq. (9), which on continum limit and after some algebra, reduces to 
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Now we modify earlier m^^(T) to m^f^{T,fi) as 

ml{T,fi) = ^nf{l + J^), (11) 

inspired by QCD perturbative calculations [2]. In our case n/ = 2, mainly due to up quarks, 
and g'^ is related to two-loop order running coupling constant, discussed earlier, but need to 
be modified to take account of finite fi. Following the work of Schneider [22] and Letessier, 
Rafelski [21], now we change T/At in Eq. (7) as 



where a is a parameter which is equal to (1.91/2.91)^ in the calculation of Schneider and 
I/tt^ in a phenomological model of Letessier and Rafelski. In our model Schneider's as{T, /x) 
works well. 

From Uq, we may obtain other thermodynamic quantities like, 

AP = P{T, fi) - P{T, 0) = ['' njfi , (13) 

Jo 

and so on. 

4 Results : 

In Fig. 1, we plotted P/T'^ Vs T /Tc for (2+l)-flavor QGP with /i = and compared with 
LGT data. We took effective number of flavors as n'^^ = 2, n'^/'l^ = .5 and n^ = 2.5 because 
of finite quark masses. Quark rest masses are m„o = Q5MeV and m^o = lQ5MeV, values 
used in LGT simulations. Only the parameter to = At/Tc is adjusted to get the best fit to 
LGT data on pressure and is equal to 0.4. In Fig. 2, the baryon density (n^/T^) is plotted 
as a function of T/Tc for different values of baryon chemical potentials and compared with 
LGT results. Surprisingly good fit is obtained without any new parameters. In Fig. 3, 
AP/T^ is plotted and again very good fit to LGT. 
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5 Conclusions : 

Using our new formulation of qQGP phenomenological model, we were able to explain LGT 
results on (2+l)-flavor QGP with just a single parameter which may be related to QCD scale 
parameter. Our formalism is thermodynamically consistent and no need of reformulation of 
SM with extra temperature dependent terms in pressure and energy density so on. Earlier, 
using this model, we explained successfully the LGT results of Bielefeld group on QGP with 
massless quarks [6]. Hence a simple model of qQGP with a single parameter, related to QCD 
scale parameter, explains all existing LGT results and may explain future results in LGT 
and relativistic heavy ion collisions. 

Of course, there are many other models which all claims to fit the LGT results, but all of 
them involve more than one parameter. Models with minimum number of parameters, which 
fits LGT results well, are SCQGP [14] and liquid model [20, 21], both with two parameters. 
All other models involve more than two parameters. Models based on QCD perturbative 
and non-perturbative calculations [12, 13] fails to fit the LGT results. Hence, it seems, 
phenomenological models based on the properties of plasma with QCD inputs explains well 
the LGT results. 
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Figure 1: Plots of PjT^ as a function of T jT^ from our model and lattice results (symbols) 
[1] for (2+l)-flavor QGP. 



Figure 2: Plots of Hb/T^ as a function of T /T^. from our model of (2+l)-flavor QGP and 
also with lattice data (symbols) [1] for /ifc = 100, 210, 330, 410 and 530 MeV from bottom 
to top. 



Figure 3: Plots of /\P/T'^ as a function of T /Tc from our model of (2+l)-flavor QGP and 
also with lattice data (symbols) [1] for /if, = 100, 210, 330, 410 and 530 MeV from bottom 
to top. 
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